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ABSTRACT 

We study the dynamical evolution of stars and gas close to the centre of the Milky Way. Any 
plausible means of forming the young stars observed at the Galactic Centre leaves behind a 
residual gas disc at ~ O.Olpc radii. We show that the combined effects of viscous accretion 
and gravitational interactions with stars do not remove the residual gas efficiently, and that a 
substantial gas disc, interior to the stellar disc, persists for > lOMyr after the stars form. Since 
no such disc is currently seen at the Galactic Centre we argue that it has been accreted by the 
super-massive black hole. This scenario offers an attractive connection between nuclear star 
formation and black hole feeding, and we suggest that the 'missing' gas may have been used 
to power Sgr A*. 
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1 INTRODUCTION 

The discovej^ of a large cluster of young stars within 0.1 pc of 
Sgr A*, the super-massive black hole (SMBH) at the centre of the 
Galaxy, poses a number of interesting questions. The physical con- 
ditions so close to the Galactic Centre (GC) are extreme, and the 
high temperatures and tidal shear prohibit the formation of stars 
by the same mechanisms that occur in the solar neighbourhood. 
However, observations show that tens of B-type stars orbit within 
O.Olpc of the SMBH, while a larger population of -100 O- 
and Wolf-Rayet (WR) s t ars exists at sligh tly larger radii, ^ 0.1 pc 
(e.g.. lGenzel et al.ll2003l: iGhez et al.ll2005l) . These stars are young 
(< lOMyr for the O- and WR stars), which suggests that they must 
have formed at or close to their current location. How stars form in 
such an extreme environment remains an open question. 

An important additional clue to the origin of the young stars 
at the GC comes from their dynamics. Many of the O- and WR 
stars are k nown to o rbit Sgr A* in a single coherent disc (e.g., 
IPaumard e t al. 20od; ILu et alj |2009| : [Gillessen et al. 200^, and 
this provides observational support for the id ea that stars can 
form in accretion discs around SMBHs (e.g ., IPaczvnskil 1 19781 : 
iKolvkhalov & Svunvaevll980l : lGoodmatj2003l) . Such discs are ex- 
pected to be gravitationally unstable at large radii, and in recent 
years a number of authors have suggested that disc fragmentation 
due to gravitational instability is responsible for the young stellar 
disc at the G C (e.g.. lLevin & BeIoborodovll2003l : lNavakshinll2006l : 
lLevinll2007l) . It is also tempting to link GC star formation to the 
broader question of SMBH accretion and growth (" King & Pringld 
[2OO6, 2007), but the relationship between nuclear star formation, 
accretion and black hole feeding remains poorly understood. 

The disc fragmentation scenario works as follows. Some sort 
of 'accretion event', such as the capture of a molecular cloud (e.g.. 
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IWardle & Yusef-Zadehll2008l) . results in gas with very low angu- 
lar momentum falling towards the SMBH. This gas dissipates en- 
ergy through shocks and radiative cooling, and forms an accretion 
disc around the SMBH. The gas at small radii accretes on to the 
SMBH, but beyond some critical radius the disc is gravitation- 
ally unstable and fragments to form stars. Numerica l simulations 
(e.g.,^ iavakshin et al. 2007; Al exander et al . 2008; Bo nnell & Ric3 
120081 - Blobbs & Navakshinll200i) ^ have shown that this scenario is 
broadly consistent with the observed properties of the GC stellar 
disc, but how the disc-star system evolves after its formation is far 
from certain. 

The fate of the re sidual gas disc is particularly interesting. An- 
alytic arguments (e.g.. lNavakshin & CuadralboOSuKing & Pringld 
I2OO7I) show that in the case of Sgr A* the accretion disc frag- 
ments into stars only at radii R > 0.03-0.05pc, which is co- 
incident with the inner edge of the observed stellar disc (e.g., 
Paumard et al. 2006; _Lu et al. 2009). Inside this radius the gas disc 
is gravitationally stable, and does not fragment into stars. Numer- 
ical simulations of gas deposition into the inner parsec find that a 
significant mass of gas, ~ IO^-IO^'Mq, is deposited into this in - 



ner region (e.g., Bonnell & Ricdl2008l : iHobbs & Navakshirill2009l) . 
and it has been argued that the residual gas disc could have been 
present for a long time, perhaps until the current epoch, supply- 
ing Sgr A* with gas at moderate rates jNavakshin & Cuadra|[200a ; 
iHobbs & Navak shin 2009). However, no such gas disc is observed 
today. An optically thick accretion disc at O.Olpc radii is com- 
pletely ruled out by observations (s ee, e.g., ICuadraetalJl2003l: 
iNavakshin. Cuadra. & Sunvaevll2004l) . and the low luminosity of 
Sgr A* itself also sets very s trict limi ts on the rate of gas accre- 
tion in any extant disc (e.g., iBower e t al. 2003). The fate of this 
'missing' gas disc remains an open question, and may provide an 
important link between nuclear star formation and SMBH accretion 
and growth. 

In this paper we model the evolution of this residual gas disc. 
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Our model is presented in Section |2] The model includes viscous 
accretion of gas and tidal torques from the stellar disc, and the stars 
interact gravitationally with one another. Our results, presented in 
Section[3] show that these processes do not operate rapidly enough 
to remove the residual gas disc in ~ lOMyr. We compare our results 
to observations and discuss possible explanations for the missing 
gas in Section|4l and summarize our conclusions in Section|5] 



2 MODEL 

Our model consists of a thin accretion disc around a SMBH, inte- 
rior to a co-planar disc of stars. The disc accretes due to viscous 
transport of angular momentum, and the stars exchange angular 
momentum with the disc through the action of tidal torques. We 
further include the interaction of the stars with one another, which 
is known to modify the ir dynamics significantly over Myr periods 
( [Alexander et al.ll2007l) . Our initial conditions are consistent with 
formation of the stars via disc fragmentation, and we use a one- 
dimensional numerical scheme to follow the non-linear evolution 
of this system forward in time. 



where Cs is the local sound speed, il = ^[GM\j^JW is the Keplerian 
angular frequency, and the dimensionless parameter a specifies the 
efficiency of angular momentum tr ansport (t hought to be due to the 
magnetorotational instability, e.g.. lBalbus & Hawle v 1991, 199^). 
We initially adopt a c onstant value a = 0.1, the fi ducial value for 
black hole discs (e.g., iKing, Pringle. & Livioll2007 ll. 

For the disc's therm al structure we use a simplified form for 
the energy equation (e.g., ICannizzd[l993h 



£7c ^ 2(6+ - Q-) 
dt c-pS 



dR 



(6) 



Here is the midplane temperature and and Q- are the in- 
stantaneous heating and cooling rates. The radial velocity ur in the 
advective term in Equation |6] is assumed to be the vertically aver- 
aged value 
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Here cp is the s pecific hea t capa city: we adopt the functional form 
for cp(T^) froml cannizzol ( Il993h . but note that the heat capacity is 
approximately constant for the low temperatures considered here. 
The rate of heating due to viscous disspation, 2+ is given by 



2.1 Disc model 

The evolution of the gas surface density, l.(R, t) in the accretion disc 
is gov erned by the equation (iLvnden-Bell & Pringle., 1974. : .Pringld 
[T98lh 
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Here t is time, R is cylindrical radius, v is the kinematic viscosity 
and Mbh = 4 x IO^Mq is the mass of the central SMBH. The first 
term on the right-hand side represents viscous evolution of the disc, 
while the second term describes the response of the disc to the tidal 
torque from the stars. A,(i?, a,) is the rate of specific angular mo- 
mentum transfer from the i\h star to the disc, and the total torque 
on the disc is found by summing over all A', stars. The torque from 
a single star on the disc is essentially t he same as that foun d in 
Type II planet migrat ion and, following iTrilling et al.l ( 1 1 9981) and 
lArmitage et alj ( l2002h . for a star of mass M, , = qiM^^^^ at radius 
(semi-major axis) a, we adopt the following form for A, 
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if R> cii 



Here 



max(/f, \R- ai\). 
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and H = cs/Cl is the disc scale-height. T his form for A, is the same 
as th at used bv lLin & Papaloizoi] ( Il986l see also lLin & Papaloizoij 
Il979i) . but modified to give a symmetric treatment inside and out- 
side the star's orbit (though in practice only the R < Oj case is 
relevant here). This transfer of angular momentum causes the stars 
to migrate at a rate 



1/2 



RA,l,dR . 



We adopt a standard lshakura & SunvaevI ( 1 19731) a- viscosity 



V = ac^H , 



(4) 



(5) 



(8) 



We adopt a 'one-zone' mo del for cooling (as used by, for example, 
Ijohnson & Gamrni3l2003h . where 



(9) 



Here, the second term allows for a smooth transition between the 
limits of optically thick and optically thin cooling, osb is the 
Stefan-Boltzmann constant, and the last term is prescribed in a 
manner than enforces a minimum disc temperature. Due to the 
strong local radiation field from the stellar cluster we expect a high 
background tem perature clos e to Sgr A*, so we adopt T^]„ = 50K 
throughout (e.g.. lLevirill2007h . The vertical optical depth t is com- 
puted as 



1 



(10) 



with the midplane density evaluated from the surface density by 
assuming that the disc is vertically isothermal: 



^/2nH ' 



(11) 



The local sound speed is essentially a means of keeping track of 
the local pressure, and includes contributions from both gas and 
radiation pressure: 



= + -; — - 



(12) 



Here K is the gas constant, fj is the mean molecular weight, and c is 
the speed of light. In the region of interest (R ~ O.Olpc) our disc re- 
mains cold, with Tc < lOOOK, and radiation pressure is negligible. 
As a result we thus fix ;u = 2.3, consistent with a disc of molecular 
gas. At these temperatures and densities the opacity is dominated 
by dust, and we use the numerical fits for k(T,p) given bv lZhu et al.l 
l l2007ll2008h . 



2.2 Stellar scattering 

In addition to feeling tidal torques from the disc, the stars also in- 
teract gravitationally with one another. Stellar interactions act to 
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increase the stellar velocity dispersion, which has the dual effect of 
thickening the stellar disc and increasing the eccentricities of the in- 
dividual stellar orbits. In practice this process is three-dimensional, 
but previous studies have shown that as long as the disc is thick 
enough to fall into the 'dispersion-dominated' regime then a sim- 
ple analytic approach i s sufficiently accurate to model the stellar 
scattering process (e.g.. [Alexander et alj2007h . For a population of 
stars with a single mass M, in a disc around a black hole, the stellar 
velocity dispersion cr, evolves as 



da, 
dt 



tr ~ CRARt„,ba-l 



(13) 



Here fj is the (two-body) relaxation time-scale, R is the mean stel- 
lar radius (semi-major axis), AR is the radial extent of the stel- 
lar disc (computed as the standard deviation of a,), and is 
the Keplerian orbital period at R. 11, = crlAR/lGM, is the usual 
Coulomb logarithnfl a nd C is a numerical scaling constant which 
[Alexander et alj ( l2007h found to be C - 2.2. 

In the absence of any gas two-body scattering causes the ve- 
locity dispersion to increase roughly as cr, oc f''*, but in the pres- 
ence of a gas disc the effects of stellar scattering act in competi- 
tion with the tidal torques described in Equations[T}H] We approxi- 
mate the effect of stellar scattering in our one-dimensional disc-star 
code by applying stochastic 'kicks' in radius to the individual stars. 
Based on the relationship between velocity dispersion and eccen- 
tricity, we define a kick size 
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where vk = ^[gMuJR is the Keplerian orbital speed at radius R. 
The stars are kicked A'k times per relaxation time t^, with the indi- 
vidual kick sizes drawn randomly from the range 



ARk ARk 
—= < oa,.k ^ — ;= 



(15) 



This form is chosen so that rate of increase of the velocity disper- 
sion is independent of the kick frequency. The value of therefore 
controls only the level of discreteness of the numerical integration. 
Test calculations indicate that the value of A'k does not strongly in- 
fluence our results, and for computational convenience we choose 
A'k = 1000. 



2.3 Initial conditions 

We choose our initial conditions to reflect the conditions imme- 
diately after a gravitationally unstable gas disc has fra gmented 
into s tars. Any plausible accretion disc model has a iToomrd 
d 19641) O parameter which de c reases with increasing radius (e.g., 
lNavakshinll2003 : lLevinll2007l ; iKing & Pringl3l2007D . and conse- 
quently disc fragmentation only happens beyond some critical ra- 
dius. We choose our initial disc surface density profile to be a 
power-law, truncated at some outer radius R^ut'- 



1,(R) oc R- 



R ^ R„, 



(16) 



and zero at larger radii. We then randomly assign the initial stellar 
radii a, in the range [Roui,2Rout]- In order to prevent Equation 1 131 
from diverging we assign a small initial velocity dispersion cr,(0) = 



' In order to avoid confusion we use H, for the Coulomb logarithm, as the 
more usual symbol (A) represents the tidal torque function here. 



20km/s [though our results are not at all sensitive to the exact value 
of o-.(O)]. 

Based on observations of the stellar discs at the GC, we model 



A'. 



100 stars of M, 



lOOMm. If these stars formed via disc 



fragmentation, the outer edge of the accretion disc should be lo- 
cated at the outermost point where the disc remains gravitationally 
stable. We therefore normalise the disc mass by demanding that the 
disc be marginally gravitationally stable at its outer radius, and set 
Q = 2 atR = i^out- In thermal equilibrium the disc temperature de- 
pends only on the orbital frequency £1 and the optical depth k, so 
fixing the value of Q in this manner uniquely determines the ini- 
tial surface density Z(i?out)- Our model therefore has only two free 
parameters: the outer disc radius R^ut, beyond which the disc is as- 
sumed to have fragmented into stars, and the power-law index p. 
Additional models were also computed in which the initial condi- 
tion was a gas 'ring', truncated inside some inner radius R,„ (see 
SectionlT2t. 



2.4 Numerical Metliod 

We solve Equations [T] and |6] using an explicit fin ite-difference 
method, on a grid that is equispaced in R'^- (e.g., IPringle et al.l 
1986), using 3437 grid cells spanning the range [0.003pc, l.Opc]. 
We use a staggered grid, with scalar quantities evaluated at zone 
centres and vectors (velocities) at zone faces. The derivative in 
the advective term (Equation |7J is evaluated as the first-order, up- 
winded value. The time-step is typically limited by the radial ve- 
locity of the gas very close to the location of the stars, and evolving 
the system on such a short time-step is unnecessary if the gas sur- 
face density here is very low. Consequently we impose a maximum 
torque per star (and therefore a maximum gas velocity in the ra- 
dial direction) of |A,| ^ OARHQ.^. We adopt zero- torque boundary 
conditions at both the inner and outer boundaries, by setting X = 
in the boundary cells. Equation[T3]is also integrated explicitly with 
the same time-step and, as discussed above, the stellar velocity dis- 
persion is modelled by giving the stars random 'kicks' A'k = 1000 
times per relaxation time. 



3 RESULTS 
3.1 Disc models 

We computed a grid of models spanning the range of plausible 
values for the two free parameters R^^t and p. The models used 
p = 1/2, 1, 3/2, 2, and R^^ = 0.02, 0.04pc, making a total of 8 in- 
dividual runs. Each of which was evolved up to a time t = lOMyr, 
which is approximately the age of the observed GC stellar disc(s). 

Fig[T] shows the evolution of the models with p = 3/2. The 
presence of the stars provides a strong tidal barrier, and the torques 
from the stars prevent viscous spreading of the gas disc to radii 
> Rou{- The back reaction of the disc on the stars is small (because 
the torque function declines very sharply with increasing star-disc 
separation; see Equation |2j, so only the inner part of the stellar 
disc is strongly affected by the gas. For the majority of the stars 
stellar scattering dominates over the star-disc torques, so the stellar 
discs spread in radius. The model with R„^t = 0.02pc evolves more 
quickly than the model with R^ui = 0.04pc, because the dynami- 
cal and viscous time-scales are shorter at smaller radii. However, 
in both cases the systems undergo only modest evolution over the 
lOMyr duration of the calculation: the stellar discs spread in radius 
by < 25%, and a significant gas disc remains out to radii > O.Olpc. 
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Evolution of model with R^jj^=0.04pc, p=3/2 
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Figure 1. Evolution of the disc-stars system: tlie upper panel shows the 
model with Rout = 0.04pc, p = Itjl; the lower panel the model with Rout = 
0.02pc, p = 3/2. In each case the lines show the disc surface density, plotted 
at ? = 0, 5 & lOMyr. The filled circles show the mean orbital radius of the 
stars, with the 'en'or bar' denoting the extent of the stellar disc (extending 
from the minimum to maximum stellar orbital radii): the top, middle and 
bottom points again represent t = 0,5 & lOMyr Stellar scattering increases 
the width of the stellar ring, while tidal and viscous torques drive accretion 
in the gas disc. However, in both cases the system only undergoes modest 
evolution: lOMyr after formation the stars have not migrated significantly 
from their initial positions, and a substantial gas disc remains. The decrease 
in surface density at very small radii is an artefact of our (zero-torque) inner 
boundary condition, and is not significant. 



These results are not substantially altered when we change 
the surface density slope p. Discs with lower values of p have 
less mass at small radii and generally evolve more quickly, but the 
overall trends seen in Fig[T]hold for all of our models. On time- 
scales ~10Myr the stellar discs do not migrate significantly from 
their initial configuration, and the combined effect of viscous accre- 
tion and tidal torques has only a modest effect on the residual gas 
disc. All of our models predict that a substantial gas disc extends 
to large radii for > lOMyr, with surface densities X > lOOgcm"' 
at i? ~ O.Olpc. These surface densities correspond to large opti- 
cal depths (t > 100), and all of our discs are extremely optically 
thick at ~ O.Olpc. By contrast, the presence of an optically thick 
gas disc at these radii is strongly ruled out by observations of the 




10^ 10^ 1 

t / yr 



Figure 2. Accretion rate as a function of time. Solid lines represent the 
models with Rout = 0.04pc; dashed lines Rout = 0.02pc. From bottom to 
top, in both cases, the curves show the models with p = 1/2, 1, 3/2 & 2, 
respectively. In all cases substantial accretion on to the SMBH persists for 
> lOMyr after the stars form. For reference, the Eddington accretion rate is 
MEdd ^O.OSMoyr-'. 

GC dCuadra et ai]|2003l : IPaumard et aDl2004) . so at first sight our 
models are not consistent with the data. 

Fig|2] shows the accretion rate on to the central SMBH as 
a function of time in the different models. We see the power- 
law behaviour characteristic of viscous accretion discs (e.g., 
iLvnden- Bell&Pringl^ll974h . modified somewhat by the stellar 
tidal torques. In all cases significant accretion on to the SMBH 
persists for tens of Myr. For reference, if we assume an accretion 
efficiency of 77 = 0.1 the Eddington rate for a 4 x 10''Mo SMBH 
is M^dd - O.OSMgyr"'. The initial accretion rates exceed O.OlMedd 
at early times, and even the lowest mass discs are still accreting at 
~ 10"''-10"''MEdd lOMyr after the outer disc fragmented into stars. 
This is in stark contrast with the current state of Sgr A*, where th e 
accretion rate is measured to be < lO^'^MEdd ( Bower et al.ll2003[) . 
Both the gas surface densities and SMBH accretion rates show our 
discs simply do not accrete quickly enough to satisfy the observa- 
tions, and we conclude that some other mechanism must remove 
the residual gas left over from the formation of the GC stellar disc. 



3.2 Additional models 

The most obvious solution to the problem of too much residual gas 
is that the initial configuration was not a disc, with gas extending all 
the way to small radii, but rather a narrow ring close to the radius 
where the stars formed. This configuration is the likely outcome if 
the disc and stars formed from an initial gas cloud with a narrow 
dispersion in its angular momentum, and allows us to place a more 
interesting limit on the post-star-formation conditions at the GC. 

In order to test this hypothesis we ran additional models where 
the initial conditions for the gas disc were a narrow 'spike' in the 
gas surface density just interior to the stellar disc. The set-up was 
as before, but now with the surface density set to zero both be- 
yond ^oui and also inside Ri„ = 0.9/?out- The surface density was 
assumed to be constant (i.e., p = 0), but we note that with such 
a narrow dynamic range in radius the surface density slope has a 
negligible effect on these models. As before we ran models with 
/?out = 0.02, 0.04pc, and this time set the viscosity parameter a = I. 
This choice of a is deliberately high, and serves as an upper limit 
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Evolution of model with narrow ring at R— 0.04pc 
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Figure 3. Evolution of the disc-stars system where the 'disc' initially con- 
sists of a narrow ring of gas just inside /?out (as desciibed in Section IT2t . 
The upper panel shows the model with R^m = 0.04pc; the lower panel the 
model with Rant = 0.02pc. As in Fig.[T] the lines show the disc surface 
density, plotted at ^ = 0, 5 & lOMyr, and the filled circles and 'error bars' 
mark the location of the stellar disc. Even with this extreme choice of initial 
conditions a significant gas disc remains after lOMyr 



to the efficiency of viscous accretion. The total initial gas masses 
in these models are 365Mq (for = 0.02pc) and 560Mq (for 
i?out = 0.04pc), which are almost negligible compared to the total 
stellar mass of 10** Mq (and imply > 95% star formation efficiency). 
This extreme choice of initial conditions is obviously highly con- 
trived, but serves to set useful limits on the evolution of the system. 

The evolution of these models is shown in Fig.|3] Despite the 
implausible initial conditions we see that the system still evolves 
too slowly. Even in these most optimistic models there is still a 
substantial gas disc present after lOMyr. The total gas mass after 
lOMyr is 342Mo in the model with R^^.^ = 0.04pc and 129Mo in 
the model with i?oui = 0.02pc. The accretion rates at the inner 
boundary are 6.9 x 10"''MQyr"' & 1.5 x 10"''Moyr"'respectively, 
and in both cases the surface density at O.Olpc exceeds lOgcm"-. 
We therefore see that even for an extreme choice of initial condi- 
tions, with an implausibly small residual gas disc and unreasonably 
efficient accretion, the combination of viscous accretion and tidal 
torques does not lead to accretion of the gas disc on a time-scale 
comparable to the age of the GC stellar discs. This suggests that 



other mechanism(s) must be responsible for removing gas from the 
region close to the SMBH at the centre of the Galaxy. 



4 DISCUSSION 

The most striking discrepancy between our model and the GC is in 
the SMBH accretion rate. There is some evidence that the accre- 
tion rate on to Sgr A* may have bee n a significant fraction of the 
Eddington rate 5-lOMyr ago (e.g., IZubovas. King. & NavakshinI 
I2OII), but the present-d ay accretion rate is very low (< lO^^MEdd, 
e.g.. Bower et ai] |2003h . We note, however, that this observed ac- 
cretion luminosity tells us only about the physical conditions very 
close to the SMBH, and does not strongly constrain the presence or 
absence of an accretion disc at larger radii. As discussed in Section 
I3.1l (see also Fig.|2]l, our models predict significantly sub-Eddington 
accretion 5-lOMyr after the stars form, ~ 10"'*-10"^MEdd- A 
SMBH disc accreting at these low rates can in principle undergo ra- 
diatively ine fficient accretion or outflow close to the Schwarzschil d 
radius (e.g., iNaravan & ^ 1 1994 iBlandford & BegelmanI 1 19991) . 
and in this case the observed luminosity of Sgr A* is not necessar- 
ily inconsistent with the presence of an accretion disc at ~ O.Olpc 
radii. Indeed, the favoured explanation for the low luminosity of 
Sgr A* is that the accretion fl ow is radiatively inefficient (e.g., 
lYuan. MarkofT. & Falckj |2002[) . Detailed study of the accretion 
flow close to the SMBH (i.e., at radii < 10^'pc) is beyond the scope 
of this paper, but we note that our models are not necessarily ruled 
out by observations of the Sgr A* accretion flow alone. 

A more stringent limit on our models comes from study of the 
disc at larger radii. An optically thick disc at radii ~ 0.01-0. Ipc 
from Sgr A* would be easily detectable, both through its intrinsic 
emission and through reprocessing of the local radiation field. At 
these radii the disc is rather cool, and the bulk of the accretion lumi- 
nosity is radiated at mid-infrared wavelengths (~ 10-100/jm). We 
can estimate the total flux emitted from the disc by using Equation 
[9] and assuming that each annulus in the disc radiates as a black- 
body. Assuming a distance to Sgr A* of 8kpc, our range of models 
gives mid-IR flux densities that range from a few mjy to several Jy 
(with larger values corresponding to more massive and more com- 
pact discs). Mid-IR observations of the inner ~ 0. Ipc are still rather 
scarce , but the recent high-resolution observations of lSchodel et"ai] 
I2OIII) show no evidence of this emission. Moreover, the presence 
of an optically thick reprocessing disc is strongly ruled out by ob- 
servations (e.g.,.Cuadra et al. 2003; Paum ard et al., 2004) . yet all of 
our models predict disc optical depths r > 10 at these radii (even 
the extreme cases discussed in Section [T2] >. Our models are there- 
fore completely inconsistent with observations of the GC. 

The basic reason for this discrepancy is easy to understand: at 
a radius ~ O.lpc, the characteristic viscous time-scale R?- Iv ~ 10'- 
lO^yr. This is an order of magnitude larger than the age of the GC 
stellar discs, so it is obvious that viscous accretion alone cannot 
remove the residual gas disc. We have shown that the combined ef- 
fects of viscous accretion and dynamical interactions with the ob- 
served stellar disc are still insufficient, and simply do not remove 
angular momentum from the gas rapidly enough. Our model is con- 
structed in the context of the disc fragmentation model for GC star 
formation, but we note that our conclusion applies to any model 
which forms stars in situ with < 100% efficiency. It is inconceiv- 
able that ~ IO^Mq of stars could have formed at t heir current loca- 
tion without creating a small resid ual gas disc (cf.. lBonnell & Ric3 
I2OO8I : iHobbs & NavakshinI I2OO9I) , but we have shown that even 
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~ IOOMq of gas cannot be accreted in a plausible time-scale. The 
key question, therefore, is what happened to this 'missing' gas? 

One potential solution would be substantially enhanced angu- 
lar momentum transport in the disc. Our models with a = i rule 
out any 'local' transport mechanism (i.e., any mechanism where 
the accretion energy is dissipated locally). Recent models suggest 
that strongly magnetised discs accrete at greatly enhance d rates 
jjohansen & LevirJl2008l : lGaburov, Johansen & Levinll201 ih . but it 
is not yet clear if this mechanism is efficient enough to drive the 
accretion required here. Alternatively, if the disc is self-gravitating 
then l ow-o rder spir al densi t y wav es can drive non-local transport 
(e.g., iBdbus & Pap aloizoj fl999l) . This is, however, unlikely to 
be significant here. The discs considered here are very thin (with 
HIR ~ 10^^), and non-local effects in self-gravitating discs gener- 
ally only become significant when the disc aspe ct ratio H/R > 0. 1 
(e.g.. lLodato & Riceil20ol ;|Durisen et &\!'200f). Moreover, the ra- 
diative coo ling time of a self-g ravitating disc at the GC is very 
short (e.g., lNavakshinll200^ : lLevinll2007h . so if the disc becomes 
self-gravitating it will fragment and form more stars, rather than 
accreting. Self-gravity can in principle explain either the formation 
of stars or enhanced rates of accretion, but not both. 

A more interesting solution is the idea that the residual disc 
could be removed via interaction with additional gas with a dif- 
ferent angular momentum vector. 'Accretion events', such as the 
capture of a gas cloud invoked in Section [T] are expected to be 
relatively frequent, and there is no observational or theoretical rea- 
son to expect successive such events to have similar geometries 
( iKinnev et aLll2000l ; [Kendall et al.ll2003l) . This picture of chaotic 
accretion has previously been suggested as a mechanism for GC 
star formation (King & Pringle 2006, 2007), and carries the addi- 
tional benefit of allowing accretion to proceed at close to the Ed- 
dington rate without sig nificantly increasing the black hole spin 
( iKing et al.ll2005l 120081) . Capture of a gas cloud with a different 
angular momentum vector can lead to cancellation of angular mo- 
mentum and very rapid accretion (e.g., Hobbs et al. 201 1), and this 
represents a plausible explanation for the missing gas disc around 
Sgr A*. Further work is needed to determine what, if any, addi- 
tional observational signatures this model would predict. However, 
chaotic accretion may explain both the presence of young stars and 
the absence of a gas disc at the GC, and provides a natural link be- 
tween the observed star formation close to Sgr A* and the larger 
question of how SMBHs accrete. 

If we accept that Sgr A* someho w accreted the 'missing' ga s 
disc, with mass ~ IO^-IO'^Mq (e.g.. iHobbs & Navakshij l2009t) . 
then the implications for the Milky Way are also interesting. Since 
this gas should have been accreted on a time-scal e much shorter 
than the age of the central starburst (^ 6Myr, e.g.. IPaimiard et al] 
[2006), then the average accretion rate on to the SMBH must have 
been a significant fraction of the Eddington rate during this puta- 
tive accretion episode. As mentioned above, this is in fact c onsistent 
with recent observational dSu. Slatver. & Finkbeinej2010l) and the- 
oretical d Zubovas et alj2011 ) evidence that Sgr A* produced AGN- 
like feedback for ~ 0. 1-lMyr at approximately the same time as the 
stellar disc formed. If this is indeed the case, the emergent chronol- 
ogy of the central parsec of the Milky Way is consistent with (1) 
deposition/capture of a large molecular cloud; (2) formation of a 
gas disc (or discs), the outer part of which rapidly fragmented into 
stars; (3) subsequent accretion of the interior remnant gas disc, 
resulting in the production of the y-ray emitting kpc-scale Fermi 
lobes. However, the detailed dynamics of these processes remain 
uncertain, and further work is needed if we are to understand this 
scenario more thoroughly. 



5 SUMMARY 

We have presented models of disc evolution following the forma- 
tion of the young stellar disc at the Galactic Centre. Our models in- 
clude viscous accretion of gas and tidal torques between the young 
stars and the disc, as well as gravitational interactions between the 
stars. We find that these processes do not remove the residual gas 
disc quickly enough, regardless of the choice of initial conditions 
and model parameters. Instead a substantial gas disc always re- 
mains, interior to the stellar disc, for > lOMyr after the stars form. 
Such a disc is completely ruled out by observations, but the fate 
of this 'missing' gas is unclear. The most plausible mechanism for 
accreting the residual gas disc is direct cancellation of angular mo- 
mentum following a subsequent accretion event, but further work 
is needed to understand the dynamics of this process in detail. Our 
results again highlights the link between Galactic Centre star for- 
mation and SMBH feeding, and we suggest that chaotic accretion 
of gas clouds provides a plausible explanation for both the observed 
stellar disc and the accretion history of Sgr A*. 
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